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ABSTRACT 


The reorder point and reorder guantity for a multi- 
echelon inventory system consisting of two levels were 
determined through the use of a mathematical model and a 
computer simulation was used to verify the results. 

The main echelon supported two lower echelon stock 
points and reordered using a continuous review inventory 
policy. The two lower echelon stock points operated under 
periodic review policies. 

The measure of effectiveness used was to minimize total 
system costs subject to a constraint on the maximum number 
of back orders per year. | 

The results from the mathematical model were used as 
input values for the simulation and measures of effective- 
ness were compared. An alternative procedure was proposed 


and simulated; and the results of the three products were 


compared. 
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fee RODUCT LON 


Despite the fact that all military supply systems and 
many large civilian corporations have multi-echelon and 
multi-location inventory systems, few multi-echelon inven- 
tory models are currently being used. The policies that 
have been used are the single echelon and single location 
inventory policies that attempt to minimize the total 
variable cost of a single location ignoring total system 
ees. Ihe Gliftficulties of Optimizing a multi-echelon and 
multi-~location inventory system are most likely due to the 
complexity of demand and the interdependency among the units 
in different echelons. 

A mathematical model for a multi-echelon and multi- 
Peete tol ventory system is developed ian Chapter If. The 
Senccm oor this medelaiseto. minimize the total annual 
variable costs subject to a maximum allowable number of 


Diemeoldememoen /ecarnefonethe vntiressystem.. This model 


cere 


consists of one first echelon location operating under a 
continuous review policy and two second echelon systems 
which operate under periodic review policies. An example 
is presented to obtain numerical results. 

In Chapter III, a computer model is developed to simu- 
late the multi-echelon, multi-location system. This simu- 


lation model was used to check the results obtained from 





the mathematical model and to demonstrate the interaction 
between echelons. Graphical plots are generated which show 
the inventory position of each entity in the model. 

In Chapter IV, comparisons are made of the simulation 
results and the results obtained from the mathematical 
model. The comparisons suggest some modifications to the 
@eerating policy. 

The system was simulated again using the modified 
operating policy and substantial improvements in the 
effectiveness of the multi-echelon system were observed. 
These results are described in Chapter V. 

Chapter VI summarizes the results of the research and 


concludes with some suggestions for additional research. 





Il. A PROBABILISTIC MATHEMATICAL MODEL FOR MULTI-ECHELON 
TMMENTORY SYSTEM 

Ae DESCRIPTION OF MODEL 

In this model there are three systems. One of these 
systems represents the highest echelon and is called the 
Main System. The other two, at separate locations, are 
the lower echelon and are called System One and System Two, 
respectively. Each system carries its own inventory, and 
receives random demands. System One and System Two are 
dependent upon the Main System but independent of each 
other. System One and System Two can only be resupplied 
by the Main System. In other words, they cannot order 
from any other suppliers. The Main System replenishes 


stocks by placing orders to external suppliers. 


Bee OBJECTIVE. OF MODEL 

The objective of this model is to minimize the total 
expected annual variable costs of three systems subject 
to a specified expected number of total back orders per 
year. In fact, this is the same as minimization of total 
yearly variable costs subject to a specified minimum level 
of customer satisfaction. 

Based on these objectives, decision variables for each 
system will be calculated to achieve optimum levels for 


the entire multi-echelon supply system. 
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C. ASSUMPTIONS OF THE PROBABILISTIC MATHEMATICAL MODEL 
1. Main System 

It 1S assumed that the Main System under considera- 
tion consists of a single installation which utilizes trans- 
action reporting. This system is governed by a (Q, xr) type 
policy with back orders. 

The assumptions in addition to the continuous review 
assumption are: 

a. The cost of operating the information processing 
system is independent of Q (reorder quantity) and r (reorder 
pewnt). 

be ane Unser cost C Of the item is a constant in- 
dependent of Q. 

c. The back-order cost is constant (II), per unit 
back ordered regardless of the length of time the back- 
orders exist. 

d. There is never more than a single order out- 
Standing. This assumption implies that when the reorder 
point is reached, there are no orders outstanding; there- 
fore, the inventory position is equal to the net inventory. 
Thus, the reorder point will be the same regardless of 
whether it is based on the inventory position or net 
inventory. 

e. Procurement lead times are independent and 
identically distributed (i.i.d) random variables with a 


gamma distribution. 


us 





f. All variables are treated as continuous. 
g. The demands are Poisson distributed with the 


mean number of demands per year a constant dan’ 


h. The reorder point, r, which is based on the 
inventory position is positive. 

With a back orders constraint, it is infeasible to 
wait until back orders exist before placing an order. 
Because of this and assumption (da) there will be no back 
orders outstanding at the reorder point. As was discussed 
in assumption (d), at the reorder point, the inventory 
position is equal to the on-hand inventory. 

For this model, any one of the three inventory 
levels on hand, net, or inventory position can be used to 
define the reorder point; and the reorder point has the 
Same value for any one of them. It should also be noted 
that to use the on-hand level it must be assumed that 
Seecer an Oraer arrives, it is sufficient to fill all back 
orders and raise the on-hand inventory level above the 
reorder point. If this ever failed to happen, the reorder 
point would never be reached again and the system would 
continue to accumulate back orders indefinitely. 

When the reorder point is thought of in terms of 
the inventory position of the system, then assumption (d) 
guarantees that the on-hand inventory will always exceed 
the reorder point when an order arrives; otherwise, it 
would not be possible to have only a single order outstanding. 
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Ze oY St emmeone 

This system, which is at the lower echelon, is 
governed by a periodic review policy. The operating 
doctrine is the most widely used type of periodic review 
which is the order up to R policy. All demands which 
occur when the system is out of stock are back ordered. 

For this periodic review system the time between 
reviews will be denoted by T, and at each review time a 
sufficient quantity is ordered to bring the inventory 
position of the system up to a level, R, regardless of the 
amount of on-hand inventory. This policy dictates that 
at review times even if the inventory position is R-l, 
only one item must be ordered to bring the inventory 
Eeotelom Up to R, igmoring the high cost of placing the 
order. Despite its appearance of being illogical, this 
assumption simplifies the formulation of the system and 
is very unlikely to occur on high demand items. When an 
item experiences zero demand in a review period, there is 
no need for order because the inventory position is 
already at R. 

The other assumptions are: 

a. The cost, J, of making a review is independent 
of the variables R and T. 

b. The unit cost, C, of the item is constant and 


independent of.the quantity ordered. 
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c. Back orders are incurred only in very small 
Manttetes. This implies that when an Order arrives, it 
1s almost always sufficient to meet any outstanding back 
orders. 

ad. The backorder cost is constant, Il, per unit 
back ordered regardless of the length of time the back 
order exists. 

e. Procurement lead times are 1.i.d random 
variables with a gamma distribution. 

£. Orders are received in the same sequence in 
which they were placed. It should be noted that for (Q, r) 
models, the two assumptions that orders were received in 
the sequence placed and that lead times are i.i.d random 
' variables could not both hold rigorously, since there 
existS a positive probability that two successive orders 
could be separated by an arbitrarily short time interval. 
In this model, orders can never be more closely spaced 
than by an interval of length T. If T is large enough, 
it is possible, provided that there is a sufficiently 
small range of variation in the lead time, that both 
assumptions hold simultaneously. 

g. The demands are Poisson distributed with 
mean, i. 


(Merle varlables ace treated as continuous. 
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3. System Two 


This system is identical to, but independent of, 


system One. 


De oNALYTICAL SOLUTION TO THE MODEL 

AS was mentioned previously, the objective of the model 
is to Minimize total system costs subject to a constraint 
on the number of back orders per year. 

The determination of various annual costs of each 
system may be made independently and placed into a common 
cost formula. Instead of doing this independently, it 1s 
better to find a main system cost expression, System One 
cost expression and System Two cost expression, and then 
to add them to each other to determine the total cost 
formulation and then use this formulation as the objective 
function of total system. The systems are actually tied 
together through the constraint on the number of back 
OmEaers. 

For the following notation, the subscript m will 
indicate that the variables being subscripted belong to 
the Main Sytem where 1 and 2 indicate System One and 
System Two, respectively. 

The costs that are of interest in this model for each 
system are the cost of placing an order, the cost of 


carrying inventory and back order costs. 
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1. Main System 
As was discussed previously, the Main System has 
a transactions reporting policy or (0, r) model. 
In a continuous review system, a period is defined 
as the length of time between the receipt of two successive 


Puscmremencs. This time period 15 a random variable. 


—— —_ 





ee 
Pacer mE = 


Because the procurement lead times are random variables, 

the number of demands for a fixed time are random variables; 
and the number of items demanded per demand are also random 
variables. The review period is also a random variable. 

In the following material, we itemize the costs. 


a. Procurement Cost 


Se = Expected number of demands per year 
oe Order Guantley 

Arn 

| = Average number of procurements per year 
m 

An = Procurement cost per cycle 

‘on 

O. An = Procurement cost per year 

m 


Peeeinventory Carrying Costs 
Because of randomness, it is possible to accumu- 
late a large number of back orders at the end of a cycle. 
To prevent this from occurring, it is advisable to provide 


a safety stock to buffer the system from excessive numbers 


ils 





of back orders. The safety stock is the expected amount 
eestock On hand when an order arrives. The actual amount 


of stock on hand when a shipment arrives is clearly random. 


Ss. = Mean value of on-hand stock when an order 
arrives 
After an order arrives, the expected on-hand inventory in- 
creases to Q+S and is reduced to a value of S on the average 
just before the next order arrives. Therefore, the average 
on-hand inventory per cycle is: 


+ 
an! ame Os 
2 ae 


f 


To write Sin in terms of the reorder point, Cn 


let us first assume that the lead time TH is fixed. bee 
— (x;r_) = r._ - xX be the net inventory 
is mM m 

at the time an order arrives, 


where xX 


is the number of units demanded in lead time Ta Then, 


= = B. [Net inventory] oe) Eom (Xi Fy) E(xiT,) Gx 
where E(x;T_) = Density function of demand in time te 
Sa pe (x 7X) £(x;T,) dx = ry a E(x;t_) dx 


= a £(x7T) ax 


ey, 





S = raTl where uw Expected lead time demand. 


Then = ~ San = ou + rm = a 
fe Inventory carrying charge 
@e— Gost Of an item 
itemal MOlaing cost/year = IC (r_ - uu. + Om 


m m 2 


Cc. Stockout Costs 
OPT: x-r <0 
m 


Were ice eleneatieles ag or nan ee 
t m ' 
m a Tee xr 20 


where n. (X >t) is the number of back orders per cycle. 


m 
ig Ney (Lan) = expected number of back orders per cycle, then 
ia Dott? “mn? hosed — i (x-r) A(x) dx 
= f xh(x) desi H(x_) where H(r) = P[X>r_] 
*m 
and h(x) = marginal distribution of leadtime demand. 


Therefore, the expected number of back orders/year = 


r 
Lf xh(x) dx - om Eis rae 


1O 
le 


and the expected cost of back orders/year = 


moe | ff och (Gx) asem cca 


18 





All the terms in the average annual variable 


cost Kn have now been found: 


i QO, HA, - 

Kn = Gu A + FC ar + we = ed ae Oo bee xh (x)dx-r_ H(r_)] 
m m m 

2. System One and System Two 


It was stated in the assumptions that System One and 
also System Two both follow a periodic review policy with an 
order up to R stockage policy. Since System One and System 
Two are identical to each other; the equations will be 
derived only for System One. 

For convenience, a period is assumed to be the time 
between the receipt of two successive orders rather than 
between the placement of two successive orders. Costs are 
described as follows: 


a. Ordering and Reviewing Costs 





J = Reviewing cost/cycle 
Pe= OTGGring Cost/cycle 
J TAY 
Ordering and reviewing cost/year = T 
i 


where T, is the period length defined in units of years. 
b. Inventory Carrying Cost 
The expected net inventory just prior to the 


arrival of an order is Ry a a AT) where U, = mean demand 


during lead time. 
The mean rate of demand remains constant over 


time and the expected demand per period must be the expected 
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amount ordered, i.e., A,T,- If the expected net inventory 
immediately after the arrival of a procurement is R)7Uy: 
it is therefore Ry Sally AST just prior to the arrival 
of a procurement. 


The expected unit years of storage incurred 


Ber period is 





aR. = me (Re ae - 4,7,)] =T,[R, - oe 

ee ee | Cae cee 1, eal ZI 
se 

and average inventory carrying cost/year = IC [Ry = ogee Vs 


Coe Seockout Costs 
First we assume the case where the procurement 
lead time is constant tT. An order placed at time t will 
arrive in the system at time t + T, and the next procurement 


will arrive in the system at time t + T+T After the 


1° 
order is piaced at time t, the inventory position of the 


system 1s R It 1S necessary to compute the expected 


1° 
number of back orders occurring between t + tT and t + T + T, - 
A back order will occur in this period under assumption c if 
and only if the demand in the time period T + T, exceeds R,- 
Assumption c also assures that after the arrival of the order 
placed at time t, there will be no remaining back orders, 

and therefore they must all occur between times t + T and 

1h ge 4E a T}. Consequently the expected number of back orders 


incurred per period is 
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cn (x-R ) if ( Xe ee T,) dx where 


eal 


ix; oT + Ty) = Demand distribution during time T + T)- 


When lead time is random with density g (tT) 
with cae and ee being lower and upper limits respectively 
and Tt, andt.,, the lead times for the orders placed at times 


Geancd t + Ty, respectively, the expected number of back 


orders incurred per period is: 


“max ‘max Bs 
ma : f Rr! (x-R,) f(x; TotT,)9(t,)9(T)) dx dt,dt, 
min min i 
“RY (x-R,) ee Ty) Ging 
i "max 
where h(x; Ty) a je (es TotT, )g(t,) dT, which is the 


Semana Gistribution during time oe Ty when lead time is 
a random variable with density function g(t.)- The average 


number of back orders incurred per year 1s; 


E, (R,,7 1? dx and the average 


ee Ore oe 


1 


beemeorder cost per year equals toll, E, (R,,T}) and 


n, (x) = Ty Ey (R) ,T.) is the expected number of back orders 


per period. 


Za 





Finally, the annual variable cost of System One 


is: 
ee oon ee lee ey (k),T,) and 





al Ty 


likewise the annual variable cost of System Two is: 


L ee 
2 : 272 
feet, Mg 





] + I, E. (R,,T,) 


where Li = J. t+ A.. 
i 7 


Be) Cbyective Function and Constraints 


The objective function of the model consists of the 


total annual variable costs of each system. Therefore; 


K = K + K, + XK which is equal to the minimization 








m i Z 
Or; 
: Q [ieee iF 
= mM _m = mms + 
Ee 5 i tC [re He = n, (4) + # 
™m m 
Nao II i eae 
eel L 2 we 2D) 
+ IC [R, a i, ] + A Ny (4) ~ T, + IC [R, Uy 
a 
+ =—— no, (r) 
T, Zz 
eee, ar) + mn, (x) + + To (x) <b 
J Qn m T) 1 T. 2 = 


where b is the specified total maximum number of back orders 


per year for the entire system. 


The problem at hand is to calculate the optimum 


values of On a Ri, ee Ro, T.5- Since the objective 


™m 


function and the constraint are non-linear functions of the 


Ze 





decision variables, we solve the problem using the Lagrange 
mMaltiplier approach. 
4. Optimum Values of Operating Variables 
After including the Lagrange multiplier in the 


formulation, the new objective function becomes: 











Minimize 
) Q a, i 
Am™ m mm - iE 
L == A +IC{(— +r - yu] + Hoe Gee 
QO, m 2 m m oO, mM tT 
A. T IT L hE 
_ Reo I 2 _ _ 2.2] 
-- Tc(R, Wy 5 ] + 7, ny (¥) + T - Ic{R, Lo 5 
I r 
2 - im = iS 1 5 
t+ a Ny x ) a Nn tt) ea ny () Tom Rs (r)) - b] 
Ez m ale Z 


where 9 is the Lagrange multiplier. It is through this 
Lagrange multiplier that the three systems are linked together 
mathematically. 

The optimum values of the unknown variables can be 
found by taking the derivatives of theobjective function 


with respect to Qn! Ra! Rye Ty. Ros To. 86; equating them to 


zero; solving the equations simultaneously; and ensuring the 


Kuhn-Tucker conditions are satisfied. 


a 


The derivatives were taken with respect to Oo 


R R 6. A different procedure was utilized to find Ty 


det oe 


and T5- 
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The derivatives: 








= ole _ ee Ler m=- 5 = 
aoe” OGm 52 ie -’ Nm 'F) (lL, - 8) = 0 
m m 
eee 
dL m 
c =—=i%1c- — (r) (TI - 6) = 0 
rn Oe On m m 
m 
ie o 
L See Pea. CRT.) UC. = 8) = OC 
Ry Ry Th i iE is ile 
= oi ogee 5 o 
L, aa Le 5 H. (Ro, T.) (I, 6) = Q 
2 2 2 
IL sae ie ie = 
se oe ee ee i) 9, (2) 
m a Z 
The Kuhn-Tucker conditions: 
Lom 2 0 Onc = 0 Q, 2 0 aan 
m 
De 2 0 Cay = 0 rn =) 8.g = 0 
m m 
Lp > 0 R,-L, an Ry a0) @< 0 
L iL 
L 2 0 Rais = Q Rea ae 
Ro 2 Ro Zz 
where g = *m n(x) + — n,(x) + ea (is) os 
oF m ee T, 2 


Therefore the following Kuhn-Tucker conditions must 


be satisfied: 
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From these equations and conditions the optimum 
values of the unknown variables can be simplified to the 
following form and the exact values can be found by 


solving these equations simultaneously: 


xn 1c oO, 
Ge) 
m eo) = 


Cn = 2r Tay + My “ly = ey ] 


ic 
. cry 
H, (Ry, T)) = cr =8) 
; Bere 

H. (Ro, T,) = Ti-8), 


The procedure to find Ty and T. and subsequently R 


and Ro PSuemrtougn 1Ceration and trial and error. This is 


a 


not unreasonable. In realistic cases other considerations 
not modelled here usually dictate the length of a period. 

In most cases the systems control a large number of different 
items and the same period length is used for each item. 

Thus, the period length is usually some convenient calendar 
or financial period, such as a month or a quarter. The 


procedure we use is shown in the following example. 


E. EXAMPLE PROBLEM 
1. Additional Assumptions 
For this example it is assumed that each system has 


Poisson arrivals independent of each other. The customers 
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demand only one item at a time. System One and System Two 
are resupplied only from the main system and the main system 
1s resupplied from outside suppliers. To be more realistic 
it is also assumed that the procurement lead-times have a 
gamma distribution with different mean and variances. 
The A's will be assumed daily demand or arrival 
rate per day rather than yearly values. 
Zee values 
Before providing the input values, some clarifications 
must be provided regarding the demand at the main system. 
When it is stated that A = 4/day, », = 3/day and x 


i 2 


These represent the mean number of demands which arrive each 


= l1/day. 


day directly at the respective systems. However, since all 
demands at Systems One and Two eventually must filter up to 
the main system, the cumulative demand at the main system 
has an expected value of dy + 5 + hin = 8. The main system 
supports not only Systems One and Two, but also has its own 
customer demands. The model assumes that the demand at the 
main system is the superposition of the direct demands at 
the main system and Systems One and Two. However, since the 
actual demands placed on the main system in the lower echelon 
are batched, the demand variability is much greater than 
would be expected by the superposition process. We will 
evaluate the seriousness of our assumption about the demand 
process at the main system with the simulation model 
described in the next chapter. 


od 





In this example the input values for each system 
are; 
a. Main System 
A. = 8/day (considering the individual customer 
arrivals and the remainder being Systems 


One and Two average demands) 


C = $50/item 
tL. = 0.23 

i = $5 

ee = $500 


b. System One 


hy = 3/day 

C = $50/item 
i = $052.3 

Ty = $5 

A, = $100 

Jy = $100 


c¢. System Two 


do = l/day 

Ge = "S507 1tem 
m= 30.23 

I, = $5 

Ay = $100 

J5 = $100 
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The total system is expected not to exceed 
b = 175 total number of back orders at the end of a year. 
3. Solution Procedure and Results 

It was assumed that for each system, procurement 
lead times are gamma distributed random variables. 

The gamma distribution is defined as: 

2 at e | with 

EAC ip 

i= a. 6 and 

Var = a. BA 

We need the lead-time demand distribution for the 
case in which lead times are gamma distributed and the 
process generating demands is Poisson. This is derived 
below. 

The Poisson distribution is 


ees est . . 
f(x) = a anGeeme demand CGdiSErIDUETON 


during lead time is: 


£ (x;T) = f(x) g(t) dt 


0 (it)* eee et /P 
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T (a+ x) dt 
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oe. 


ae 


a(l- Wine 
ee ee ao eee = AaB 


1/8 +A 


1/8 
a (l-+7ery) 
a 
1/8+A 


5 = AaB(1+AB) 


c 


a Negative Binomial distribution 
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The parameters a and 8 for each system are: 


Cn = 12.8 ay = 5.43 a, = 2.5 


i Bek 25 By 


BW 
it 


The main system has a Negative Binomial distributed 


lead time demand with parameters: 


peo eo) 5.125)" = 9320 


var. = 8° (12.8) (3.125) [1+8-(3.125)] = 8320 


system One has the parameters: 
Nee Boo. 44.6) =2 74.954 


eer Ss (9543) (4256) (143 (426),) = 1209. 0232 


i 
System Two has the parameters: 
Hy = i225). 6 = 15 

var, = Hz) oma lL + (6) | — 105 

Because the negative binomial is computationally 
intractable, we use the normal approximation for the calcu- 
lation of lead time demand probabilities. The normal 
distributions are assumed to have the same mean and variance 
as the negative binomial distributions they replace. 


To solve this example problem, it is necessary to 


make an initial estimate for the Lagrange multiplier § 
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which must be less than or equal to zero. After the initial 
estimate if the total number of back orders per year exceeds 
b, then the absolute value of 6 should be increased gradually 
until the number of back orders converges to b. 

For this problem ®@= -1.5 works very well. For the 


main system the formulas are: 


oO, Yea, l4, * Ny te) C8 Ic 
on EC 

H (xr) = = C= CP SC[ XX er] 

m (Il Oo) A, m 


2. is calculated to be 550 and 


rf Be s50 Or 3) 50 
Hote) = P [X2z,,] = Sec—eygonp = 0-3332455216 


From the inverse standard normal distribution TI-59 
calculator program (Appendix G), the reorder level at the 


main system is found to be: 


Mo om GO 725053353395) + . = 35922706827 
m m m 


‘mm, 
S 
m 





A,r) = (Hon? Ht) BOs B ( 


where @ is the functional value of standard normal distribu- 
elem at 0.4305333395 


ee) = 20.08139334. Using this value in the Q., 


formula yields 


Bi 





= 72.292 : 
on San 2 OSS ISSA (SPS) en: ag 


(ORZey 5.0 
which is not equal to the first estimated value. If on is 


taken 569 then: 


Peso. co) 50 
H (x) = (errs eae = 0.3447576396 


a. COs sI9071 201) = i 356.4008941 


n (©) = 21.05458646 


Using this value in the Qn formula yields 


Q = ¥2.2920 [500 + 21.05438646 (5 + 1.5) = 568.69 


Oma 


which is very close to the initial o value. A summary of 
the results for the main system follows: 


= = 569 


> 36.400891 


mm 


n(x) = 21.05438646/period 
n(x) =108.047115/year 


Kn = $6796.50/year 


To solve for the optimum values for Systems One and 
Two a different procedure is followed. The total annual cost 
of Systems One and Two is a convex function of the period length 


T. For different values of T there are different total annual 
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cost values. The minimum of these values is the optimum 
total annual variable cost and the corresponding T and R 
values are the optimum operating values. A TI-59 program 
waS written to perform the line search for the best value 
of T. The program is found in Appendix G. User information 
and the features of the program are also in the same appendix. 
This program evaluates the R value, safety stock, back 
orders per period, back-orders per year, annual reviewing 
and ordering cost, annual holding cost, annual back order 
cost, and finally, the total annual cost. 

Using this program the optimum values for Systems 


One and Two are: 


Ry = 306.809389 

Ty = 2.39 months (72.6958 days) 

34 = 13.78788901 

n(x) = 8.670289086/period 

fiy(x) = 43.53283223 

K, = $2634.41/year 

R, = 134.9807003 

T., =A Or montis sles. 7t5e5535 days) 
Si = 0 

No(2) = 8.13480099/period 
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fee) 23.98467122/year 


Ko 


$1377.56/year 


The total back orders per year are: 


n(x) + 1, (2) + no (x) = 108.047 + 43.532 + 23.984 


l= 56 


which is the same as b = 175.56. 
The total annual variable cost is: 


K = Kn + Ky ts Ke =o 90, 08+) coseeol + ilo). 56 


A 
il 


JOS Cb.4 7. 


Bis: 





ITI. COMPUTER SIMULATION FOR MULTI-ECHELON MULTI-LOCATION 
SINGLE ITEM INVENTORY SYSTEM 

To check the analytical results, a computer simulation 
was written which uses the same operating assumptions and 
input parameters that were made for the analytical model. 
This model, however, simulates the real world more accu- 
rately since some of the simplifying assumptions required 
to obtain analytical results were not necessary in the 
Simulation. Also, the demands placed at the main system 
from the lower echelon systems were batched as in the real 


Wwe id . 


Aw DESCRIPTION OF THE SIMULATION MODEL 

As in the analytical case, there are three systems in 
the simulation model; the main system, System One and System 
Two. The flow charts of this program are in Appendix A. 

Pec System 

This system uses a continuous review policy. When 

the stock level reaches the reorder point, it orders the 
amount Q. The decision variable for reordering is the in- 
ventory position. When an order is placed, the inventory 
position increases by an amount of Q. If the new inventory 
position is less than the reorder point, the system places 


an additional order for another amount of size Q. Then the 
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inventory position increases with one more Q. This con- 
tinues until the inventory position exceeds r. The order 
policy is thus (nQ,r). N is the smallest integer that will 
make the inventory position higher than the reorder point. 

There are independent customer arrivals to the 
Main system and also group demands that are placed by 
System One and System Two when those systems replenished 
their stocks at the end of their periods. 

The number of units demanded per requisition is 
one. The number of units demanded for resupply to the lower 
echelon systems is random depending on the demands at the 
L? Ro: The 
program is capable of allowing geometrically distributed 


lower echelon systems and the parameters R 


quantities demanded per requisition. Back orders at the main 
System are satisfied by filling the back orders to individual 
customers first followed by filling any back orders due to 
System One and System Two on a first-come, first-served 

basis. 

When a demand occurs from the lower echelon systems 
for which there is not sufficient on-hand inventory at the 
Main system, the maximum amount is filled and the rest is 
put into the back order queue. As soon as a shipment 
arrives at the main system, the back orders are filled. 

The times between the customer arrivals are independent 
of each other and exponentially distributed. The lead 
times are also independent of each other and gamma distributed. 
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2. System One and System Two 

Both systems operate identically. The only differ- 
ence might be in the values of the system parameters. They 
have the same periodic review policy which is: order up to 
R at each review time. The program is also capable of making 
a decision at each review time regarding the placement of 
an order if the inventory position is less than or equal to 
Beenresnold value, is an (r,R) policy. If r is taken to be 
R-1, then the system orders up to R at each review time even 
if there is only one demand in a period. 

The times between customer arrivals are independent 
of each other and exponentially distributed. The number of 
units demanded per requisition is one. The lead times are 


independent of each other and gamma distributed. 


feeeiewetO USE THE PROGRAM 
Two Simulations are given in Appendix C and Appendix D. 
The programs differ primarily in the type of output that 
is generated. The program given in Appendix C produces a 
Versatec plot output showing the inventory position for 
the main system and the lower echelon systems. A simulated 
period of up to 4 or 5 years can be run with a Class K job. 
On the Versatec output for each system, there are two 
plots: one plot has a small triangle at each point repre- 
senting the inventory poSition. The second plot represents 


the net inventory. The first figure shows the main system 
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values, the second figure shows the System One values and 
the, third figure shows System Two values all plotted versus 
forme 

The second program can be used to simulate operations 
over arbitrarily long periods of time. There are no dimen- 
Sion restrictions. This program gives the net inventory 
at the end of each period for all the systems; total demand 
in the periods of both System One and System Two; demand 
during a lead time for the main system; average on-hand 
inventory; and average number of items short per day for 
each system. 

The input required by the programs is described by the 
variable definition list given in Appendix H. The programs 


use the IMSL subroutines for random number generation. 


e-) RUNS 
IMeccortaing Conditions 

There were two different starting conditions entered 
for the main eee In the first case the inventory posi- 
tion and net inventory were set initially to equal the order 
quantity. In the second case, the inventory position and 
the net inventory were set equal to the reorder point plus 
the order quantity. For Systems One and Two the inventory 


positions and net inventories were set initially at R,. 


Se 





De Results of Runs 


In all runs the length of time was 10 years. The 


average yearly results were obtained by dividing the results 


for the 10 year period by 10. 


Using the policy parameters determined by the mathe- 


matical model and the first set of starting conditions the 


following simulated results were obtained. 


ae 


Main System 

The number of orders in 10 years = 5l. 

The number of back orders in 10 years = 3026. 
The average on-hand inventory over 10 years = 
Ble O26. 

The average safety stock when an order arrives = 
weno 4 un. tS 


The average annual costs are: 


Ordering Costs: S27 550800 
Holding costs: 3,599.80 
Back Order Costs: i oes . 00 


TLObaweCost 57 7662.00 
System One 
Number of back orders in 10 years = 587. 
The average on-hand inventory over 10 years = 


ita | ees Un te sy 
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The average annual costs are: 


Ordering and Reviewing Costs: $1,004.18 


Holding Costs: eo One 
Back Orders Costs: 293.50 
Total Costs ee, 614.22 


c. System Two 
Number of back orders in 10 years = 287. 
The average on-hand inventory over 10 years = 
De oom unl tS 
The average annual costs are: 


Ordering and Reviewing Costs: S$ 589.68 


Holding Costs: 5G2 2.01 
Back Orders Cost: 143.50 
Total Costs Si 39529 


(1) Total results 
The total number of back orders for the 
whole system per year is 302.6 + 58.7 + 28.7 = 390. 
The total annual variable cost for the 
TIOUeovyceen 1S 57,662.80 + 52,614.22 + $1,395.79 = $11,672.81. 
The model was run again with the second set 
of starting conditions and the same set of policy parameter 
values. The Simulation results are summarized below: 
(a) Main System 
The number of orders in 10 years = 50. 
The number of back orders in 10 years = 


EGO or 
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(b) 


(c) 


The average on-hand inventory over 

10 years = 346.937. 

The average safety stock when an order 
arrives = 77.46. 


The average annual costs are; 


Ordering Costs: S290 6,010 
Holding Costs: Bro cine 
Back Order Costs: S32eoe 
Total Costs: SB? O22 se 
system One 
Number of back orders in 10 years = 582. 


The average on-hand inventory over 
10 years = 118.696 units. 
The average annual costs are; 


Ordering and Reviewing Costs: $1,004.18 


Holding sGoscs: Ne 6s00 
Back Orders Costs: 29 aoe 
LOtcalweOstes <= : - SZ), oe Oss 


System Two 

Number of back orders in 10 years: 314. 
The average on-hand inventory over 10 years 
57.401 units. 


The average annual costs are; 
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Ordering and Reviewing Costs: $ 589.68 


Holding Costs: 660.11 
Back Orders Costs: ie 0 0 
Total Costs: e406. 79 


Ca Total Results 

The total number of back orders for whole system 
per year is 166.5 + 58.2 + 31.4 = 256. 

The total annual variable cost for whole system 
omy eoe2.2ao + 22,600.18 + $1,406.79 = $11,389.25. 

3. Comparison of Both Starting Conditions 
It 1S obvious that if both results are compared there is 

a considerable decrement on annual back orders but not a great 
difference in annual variable costs. 


B90-—256 


390 = 0.3436 


ogee = 11389325 


Layer oul = 0.02429 
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IV. COMPARISON OF ANALYTICAL AND COMPUTER SIMULATION RESULTS 


By comparing the analytical results with the simulation 
results, it is possible to evaluate the analytical model. 
If the measures of effectiveness predicted by the analytical 
model are reasonably close to those generated by the simula- 
Pion: meme 1S Support for the analytical results. in the 
table below we compare the main measures of effectiveness: 


back orders and costs. 


Total Back Orders Total Variable 


Per Year Costs Per Year 
Analytical result i's. 6 SLOFS08 .17 
Simulation result 390.00 Sul owize 6 |) ote 
(For both starting 
conditions) 25620) Slivs3e9.25 


This table shows that the average yearly back orders are 
considerably higher than what is estimated by the mathematical 
solution. In fact, the total number of back orders generated 
in the simulation does not even meet the constraint. Since 
the solutions produced by the simulation model, especially 
with respect to the number of back orders do not agree with 
the mathematical model, this suggests that some of the 
assumptions made in deriving the mathematical results are not 
reasonable. The major differences between the simulation 


results and the analytical results are found in the measures 
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of effectiveness at the main system. The results for the 
lower echelon track reasonably closely. 

MMemCOstsS are also a little higher than what is 
expected but much closer percentage wise. 

AS an explanation for the large differences observed 
in the number of back orders at the main system with the 
Simulation and analytical models, let us reexamine the 
assumptions we made that affect back orders. Back orders 
results from inadequate amounts of safety stock to protect 
the inventory system against excessively large numbers of 
demands in a lead time. The safety stock is manipulated 
iaeeconerol Of the reorder point r. After the reorder point 
is hit and an order placed, the system is totally at the 
mercy of the demands that occur during the lead time. If 
the variance of lead time demand is underestimated, then 
large stockouts will occur, even if the mean lead time 
demand is estimated accurately. The lead time demand is 
affected by two random quantities: (1) the demand distri- 
bution and (2) the lead time distribution. 

For purposes of making the mathematics tractable, we 
assumed in our analytical model that demands at the main 
system flowed in at a smooth continuous rate A which was 
taken to be the sum of the rates of demand incurred directly 
at the main system a and those which occurred at the lower 


echelon systems Ay and o> The lead times were assumed to 
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be normally distributed with mean equal to At and variance 
- \.(1- Ag) where t is the mean vate Sp veiw Ieee! eaves 
Because of the relatively high demand rates used in the 
example runs, the normal assumption should be justified by 
the Central Limit Theorem. However, as is illustrated by 
the Versatec plot output (See Appendix F) the demands at the 
main system are far from smooth. What happens is the demands 
which arrive directly at the main system cause the inventory 
position to drop off smoothly. However, when the replenish- 
ment orders from the lower echelon are received, large drops 
occur in the inventory position of the main system. Recall 
that the lower echelon systems order in batches once each 
period from the main system. If the main system simply tries 
to average out demands (as assumed by the analytical model), 
it will sometimes have very large amounts of excess stocks 
when shipments arrive and sometimes very large numbers of 
back orders. The high variance in lead time demand caused 
by the irregular demand actually seen at the main system 
causes the problem. Since all demands eventually flow 
through the main system and this is assumed by our analytical 
model, the problem cannot lie in the value used for the 
mean lead time demand. 
Let us explore further what happens when reorders are 
triggered at the main system. Demands directly at the main 


system eat away smoothly at the inventory position. Then a 
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very large quantity, say X, is demanded by one of the lower 
echelon systems. There is a very good chance that the 
large order placed by the lower echelon system will trigger 
a reorder by the main system. However, if the demand causes 
a large overshoot of the reorder point, the main system may 
awemmuecn less stock to live off of until the shipment 
arrives. For example, suppose the inventory position is 
ener, ee reorder point is r = 300, and System One places 
it's resupply order for 200 units. An order will be placed 
by the main system but instead of having 300 units to keep 
it going until the order is received, it will have only 
iwmenrts., Lt 1s clear that if 300 is the amount of stock 
needed to provide reasonable protection against demands in 
a lead time, large numbers of back orders would be expected. 
Effectively, in the example above, the reorder level was not 
r = 300 but r = 147 and the safety stock negative. The 
impact of this surge in demand caused by the batching of 
demands received from the lower echelon systems is to reduce 
the "effective" reorder point from the value r to a value 
Maeeeesccause the actual demand distribution witnessed 5, 
the main system is difficult to describe mathematically, 
the actual value of r' cannot be determined analytically. 
However, it is clearly less than r and may be much more so. 
In the next chapter we describe a policy modification 


for operation of the main system that was suggested by the 
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observations above. The modification attempts to allow 
the main system to anticipate the surge of demands that 


will be received by the lower echelon systems. 
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V. AN ALTERNATIVE SOLUTION 


A. DESCRIPTION OF ALTERNATIVE PROCEDURE 

Due to the reasons mentioned in Chapter IV, the number 
of back orders found by simulation were much higher than 
the number of back orders predicted by the mathematical 
model. The question is how to run the multi-echelon system 
so that the large number of back orders seen earlier can 
be reduced. 

Obviously, what needs to be done is to reduce the impact 
of the very large demands that occur when the lower echelon 
systems place their resupply orders. With modern day 
communication and data systems, it would be feasible to 
allow the main system to "See" every demand that occurs 
anywhere in the system. If the main system is given the 
visibility, it will be able to take action to get the stock 
on its shelves in anticipation of the large demands from 
the lower echelon systems. The main system can do this if 
it uses as its reorder point the pseudo inventory position 
which is like the inventory position except that it decreases 
only when direct customer demands are encountered at any 
of the three systems. The pseudo inventory position is un- 
affected by the batch replenishment demands placed by the 


lower echelon systems. For example, if a customer requests 
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a unit from System One, the pseudo inventory position at the 
main system decreases by one. However, if System One places 
an order for 200 units at the main system, the pseudo inven- 
tory position does not change. 

This new policy is referred to in this thesis as the 
"early warning policy." Clearly, the pseudo inventory posi- 
tion will always reach the reorder point before the inventory 
position. Therefore, orders will always be placed earlier 
and consequently, the number of back orders should decrease. 
The price paid will be in terms of extra holding costs. The 
results of a simulation using the early warning policy should 
be more nearly like those of the mathematical model since, 
in effect, the mathematical model makes the assumption of 
early warning. Implicitly, assumption of demand which is 
the superposition of the direct demands occurring at the 
main system and Systems One and Two is equivalent to the 
"early warning assumptions." 

In the next section, results are given of a simulation 
of the multi-echelon system with early warning. The flow- 
chart of the simulation model is given in Appendix B. The 


actual FORTRAN computer program is given in Appendix E. 


B. COMPUTER SIMULATION RESULTS FOR EARLY WARNING POLICY 
The same starting conditions and input parameters used 
in the previous simulation run were utilized here. The 


results are summarized as follows: 
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c.. 


yearly back orders and total variable system costs 


1. Main System Results 


The number of orders in 10 years = 50. 


Total number of back orders in 10 years = 126. 
The average on-hand inventory = 488.093. 
mwetcage Salety stock = 210.96. 
The costs are: 
Ordering costs = $2,500.00 
Holding costs = 5,613.00 
Back wOrgder Costs = 63.00 
MOEal costs S38 176 07 
2. System One Results 
Total number of back orders is 10 years = 665. 
Biceavelege On=Nand inventory over LO years = Il7. 346, 
The costs are: 
Ordering and reviewing costs = $1,006.18 
Holding Costs =e 5 es 
Bacrorder Costs = 33:2 256 
Total costs S2,691.9on 
3. System Two Results 
Total number of back orders in 10 years = 312. 
The average on-hand inventory over 10 years = 58.552. 
The costs are: 
Ordering and reviewing costs = $ 589.68 
Holding costs 3 67355 
Bacmmoncdens COSTS = 156.00 
Total costs Sl, 4930s 


COMPARISON OF EARLY WARNING POLICY SIMULATION RESULTS 


WITH ANALYTICAL AND FIRST SIMULATION RESULTS 


The comparison will be done with respect to 


Sidering all three system entities). 


ou 


total 


(con- 








Entire system 


Number of Back Total. 
Orders Per Year Cost 
Analytical Result 175256 meoreors 08.47 
Pesto imulation Results 256 ieee > 
Early Warning Simulation 
Results Troe Eee 2 Ok 


Since the effects of implementing the early warning policy 
will be observed primarily at the main system, we also 


Peeduce the results obtained for the main system individually. 


Main system 


Numbers of Back Total 
Orders Per Year Cost 
Analytical Result 108.05 S65 96450 
First Simulation Result 166.50 7,322.25 
Early Warning Simulation 
Results eZee 8 le Ow 


This table shows the differences better than the first one. 


The number of back orders decreases 88.3 percent, simultaneously 


as the costs increase about 10.44 percent. The reason for the 
higher cost is becausethe increment in safety stock increases 


the carrying costs. 
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VI. CONCLUSIONS 


The differences between the resultsof the mathematical 
and the simulation model can be explained largely as a result 
of the assumptions made about the demand process. In the 
main system it was assumed that the demand was smooth (the 
superposition of three Poisson processes), but in the simu- 
lation model the actual demand at the main system was as it 
would be in actual practice. There were batches of demand 
placed by System One and System Two in addition to the indi- 
vidual customer demands directly at the main system. These 
demand batches in fact increased the variance in lead time 
demand beyond that modelled. This explains why many more 
back orders were generated in the simulation model than 
what was predicted by the mathematical model. 

The simulation model developed in this thesis is useful 
for making comparisons and examining the effects of policy 
changes or parameter changes. Moreover, it is one of the 
best ways to check the reasonableness of all of the simpli- 
fying assumptions made in order to obtain analytical solutions. 

The mathematical results described in this thesis do 
not adequately determine the reorder point. The predicted 
number of stockouts is much less than the simulated numbers. 


As explained earlier, this is probably due to the assumptions 
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made in the model about the variance of lead time demand. 
Further study is needed to determine what could be done 

in the analytical model to better approximate what happens 
timeactual practice. 

The early warning policy discussed in this thesis did 
provide for great reductions in the number of stockouts 
system wide. Since stockouts are probably the most important 
consideration in military supply systems, the early warning 
policy 1S recommended, even though the holding costs are 
larger. Additional study is required to see if the reorder 
level in the early warning policy can be reduced substantially 
from the value determined by the mathematical model. Pre- 
liminary evidence is that the reorder level can be reduced 
significantly (25 percent or so) without generating excessively 
many back orders if the early warning policy is used. 

Our simulation models allow us to view the effect of 
changes but they cannot be used to optimize the values of 
the policy parameters. For that objective, additional work 
in the mathematical modelling area 1s required. 

fimeenis thesis we have tried to model a multi-echelon 
inventory system analytically, by linking together the 
individual echelons and locations through a single objective 
function and a constraint on backorders system wide. We, 
knowingly, were making various simplifying assumptions to 


facilitate the derivation of solutions. As reported above, 
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the resulting solution for the reorder level at the main 
system led to many more back orders than predicted. The 


other solutions; R and T, appear to be satis- 


le Ros Ge The 2 


RACTOLY . 

In order to model adequately the multi-echelon system, 
1t will probably be necessary to build into the determina- 
tion of the reorder lead at the main system the values of 
the parameters Ri and Ro at the lower echelon systems. 

This will be the only way to accurately describe the actual 
demand process that is observed at the main system. We 


recommend future work in this area. 
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APPENDIX A 


PLOWCHARTS OF FIRST SIMULATION MODEL 


Subroutine One = Ship arrivals to System One 
Subroutine Two = Ship arrivals to System Two 


Subroutine Three = Ship arrivals to Main System 


Subroutine Four Periodic review of System One 


Subroutine Five Periodic review of System Two 
Subroutine Six = Shipment arrival to System One 
Subroutine Seven = Shipment arrival to System Two 


Subroutine Eight = Shipment arrival to Main System 
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APPENDIX B&B 


Flowcharts of Early Warning Simulation Model 


Subroutine 
Subroutine 
Subroutine 
Subroutine 
Subroutine 
Subroutine 
Subroutine 


Subroutine 


One = Ship arrivals to System One 
Two = Ship arrivals to System Two 
Three = Ship arrivals to the Main System 


Four = Periodic review of System One 
Five = Periodic review of System Two 
Six = Shipment arrival to System One 
Seven = Shipment arrival to System Two 


Elght = Shipment arrival to the Main System 
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APPENDIX D 
Simulation Program; No Versatec Plotter Output 
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APPENDIX F 
Versatec Output of the Program in Apvendix C 
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APPENDIX G 


iti-=sJeecaleulator Program Eor Analytical Solutions for 


Periodic Review Systems 


This program was written for probabilistic periodic 


review inventory models having gamma distributed lead times 


and Poisson arrivals and uses. Normal distribution for 


demand during lead time instead of Negative Binomial. 


Input Requirements 


The following variables should be stored in the regis- 


ters shown before the variables. 
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also requires the TI-59 applied statistics 
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APPENDIX H 


Variable Definitions for Simulation Programs 


il 


Net inventory 

Inventory position 

Second inventory position for early warning system 
Order quantity 

Reorder level 

Number of items demanded by systems 

Number of items demanded per demand by ships 

Ship queue 


Lower echelon's group demand queue in the main 
system 


Amount of demand for each demand waiting in the 
queue 


Amount of shipment arrived 

Index. If it 1S equal to l, that means that the 
demand waiting in the main system queue to be 
filled belongs to System One. 

Total waiting time 

Total average on-hand inventory 

Length of a period 

This indicates the subroutines 

This indicates time of subroutines scheduled 

Net inventory variable for Versatec plotter 
Inventory poSition variable for Versatec plotter 


Time for Versatec plotter 
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WS 


WZ 


IK 


WW 


SS 


IQQ 


ML 


KA 


py 


Work space for geometric random variable 
Work space for geometric random variable 
Work space for geometric random variable 
Exponential random number 

Geometric random number 

inerement 


Indicates the number of ships waiting in the ship 
queue at System One 


Indicates the number of ships waiting in the ship 
queue at System Two 


Indicates the number of ships waiting in the ship 
queue at the main system 


Indicates the number of demand batches waiting 
in the group demand queue at the main system 


Time indicator 


Counter for ship arrivals to indicate the number 
of items demanded per period 


Multiplier for the number of batches of demand to 
be ordered from outside supplier 


Indicates the last change on net inventory of the 
main system in order to get the number of items 
demanded in a lead time 

Time of last change 

Scale parameter for a gamma distribution 

Shape parameter for a gamma distribution 


Ship arrival rate per day 


Probability of success for geometric distribution 
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Healey ,.G., and Whitin, T.M., Analysis of Inventory Systems, 
Prentice-Hall, 1963. 
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